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Abstract 

We study the rectification process of interacting quantum particles in a periodic potential exposed 
to the action of an external ac driving. The breaking of spatio-temporal symmetries leads to 
directed motion already in the absence of interactions. A hallmark of quantum ratcheting is the 
appearance of resonant enhancement of the current (Europhys. Lett. 79 (2007) 10007 and Phys. 
Rev. A 75 (2007) 063424). Here we study the fate of these resonances within a Gross-Pitaevskii 
equation which describes a mean field interaction between many particles. We find, that the 
resonance is i) not destroyed by interactions, ii) shifting its location with increasing interaction 
strength. We trace the Floquet states of the linear equations into the nonlinear domain, and 
show that the resonance gives rise to an instability and thus to the appearance of new nonlinear 
Floquet states, whose transport properties differ strongly as compared to the case of noninteracting 
particles. 

PACS numbers: 05.45.Mt, 05.60.-k, 32.80.Pj, 05.45.-a 
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I. INTRODUCTION 



The breaking of space-time symmetries, and their role in the generation of directed trans- 
port in single particle Hamiltonian ratchets, have been extensively studied in the classical 
pl I2I, [J and quantum regimes 0, 0, [(j, [9, Q. Experiments with thermal cold atoms loaded 
on optical lattices |tJ demonstrated the fruitfulness and correctness of the theoretical pre- 
dictions. Importantly, the latest studies show a resonant enhancement of the current in the 
quantum regime, due to resonances between Floquet states [5|, |6| . Real experiments involve 
many atoms, and interaction between them may be tuned, but will always be left at least 
at some residual nonzero level. Therefore, the impact of interactions on quantum ratchets 
has to be addressed. 

In this paper, we study, using the mean-field approach, the generation of directed transport 
of interacting quantum particles in a periodic potential under the action of a two harmonic 
driving. With this approach we mimic the motion of cold atoms in an optical lattice under 
the presence of an external force 7| , but at much lower temperatures, when a Bose-Einstein 
condensate may form To this end, we investigate the continuation of Floquet states 
of the corresponding linear system into the nonlinear domain. We show that a resonant 
enhancement of the current in the nonlinear regime takes place, which results from the 
resonant interaction between Nonlinear Floquet states. We derive an analytical expression 
for the evolution of quasienergies in the nonlinear regime. Finally we show the relation 
between the transport properties of nonlinear Floquet states and the asymptotic current of 
an initial state with zero momentum. 



II. MODEL 



Experimental realizations of ratchets with cold atoms may tune the temperatures from mK 
down to /zK, such that a Bose-Einstein condensate may form due to interactions between 
particles [8|. The corresponding general equation to be studied is then given by the one- 
dimensional Gross-Pitaevskii equation (see e.g. Ill ]) 
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where a s is the s-wave scattering length, M is the atomic mass, = ir/d is the optical 
lattice wave number with optical step d, V is the periodic potential depth, and e(r) is a 



periodic driving force. The wave function is normalized to the total number of atoms in the 
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condensate and we define uq as the average uniform atomic density 
Introducing the dimensionless variables x = 2k L X, t = r/t s , ip = ^ / \/no, and defining 
1/yU = M I %fik\t s ] we transform the system (CD) to the dimensionless equation [3] 



where the dimensionless one-particle Hamiltonian is 

1 



7:p 2 + v cos(x) - xE(t), 



(2) 



(3) 



with p = —i^jlj: and rescaled parameters v = /2 2 MV / '&h 2 k 2 L = 1, E(t) = jj 2 Me(t)/8h 2 k 3 L 
and g = [i 2 C with C = imoa s /k\ The dimensionless ac field E(t + T) = E(t). 
As in the linear limit 0, [f3], we consider E(t) = E\ cos[c<j(t — 1 )] + E 2 cos[2u;(t — to) + &] with 
to as initial time. 

By using the gauge transformation, — > exp[^xA(t)] where A(t) = — E\ sin[a;(t — 
t )]/uj — E 2 sm[2u(t — to) + 0]/2lo is the vector potential [5|, [6J; we transform the original 
one-particle Hamiltonian to 



H = -[p-A(t)] 2 + cos(x). 



(4) 



A. Linear regime 



Consider the Schrodinger equation, the linear limit of Eqs. 



Here, we use a tilde to 



denote wavefunctions, quasienergies and other relevant parameters for the linear regime. 
It was shown in 5|, |6| that, for the appearance of a dc-current in the quantum regime, two 
symmetries need to be broken. These symmetries are defined in the classical limit as follows 
PJ. If E{t) is shift symmetric E{t) = —E{t + T/2), then the Hamiltonian ([3]) is invariant 
under the transformation 



S a : (x,p,t) -> (-x, -p,t + T/2). 



(5) 



Likewise if E(t) possesses the symmetry E(t) = E(—t), then ([3]) is invariant under the 
transformation 

Sb ■ (x,p,t) -> (x, -p, -t). (6) 



The Hamiltonian Eq.([3]) is a periodic function of time. Then the solutions, \ip(t + to)) — 
U(t, t )\ip(to)), can be characterized by the dynamics of the eigenfunctions of U(T, t ) which 
satisfy the Floquet theorem: 

\Mt)} = e-^lMt)}, \4> a (t + T)> = \} a {t)). (7) 

The quasienergies e a (— tt < e a < tt) and the Floquet eigenstates can be obtained as solutions 
of the eigenvalue problem of the Floquet operator 

U(T,t )$ a (t )) = e-^\Mt )). (8) 

Due to the discrete translational invariance of Eq.(T4]) and Bloch's theorem all Floquet states 
are characterized by a quasimomentum k with \ij) a (x + 2tt)) = e lhK \ip a (x)) . 
We choose k = which corresponds to initial states where atoms equally populate all (or 
many) wells of the spatial potential. This allows us to use periodic boundary conditions for 
Eq.(j2J), with spatial period L = 2it, so that the wave function can be expanded in the plane 
wave eigenbasis of the momentum operator p, \n) = -^=e mx , viz. 

N 

\m) = J2 tn(t)\n). (9) 
n=-N 

Thus, the Floquet operator is obtained by solving Eqs.(j2]HD in the linear limit. In the 
computations we neglect the contribution originating from A(t) 2 , since it only yields a global 
phase factor. Details of the numerical method are given in [6] . A study of a related problem 



with nonzero k has been published in which shows that the essential features of avoided 
crossings survive. Once the Floquet evolution operator is computed, one finds that the 
symmetries of the classical equations of motion are reflected by corresponding symmetries 
of the Floquet operator. If the Hamiltonian is invariant under the shift symmetry S a 
then the Floquet operator possesses the property U(T,t ) = U* T (T/2,t ) U (T/2,t ) 
Here performs a transposition along the codiagonal of U . With our driving function E(t), 
S a is always violated. If the Hamiltonian is invariant under the time reversal symmetry S& 
(EJ, then the Floquet operator has the property U(T,t ) = U{T,t )* flflQ. 
To estimate the net transport, it is necessary to compute the asymptotic current. It is 
obtained using the expression J(t ) = J2a(p)a\C a (to)\ 2 [6], where (p) a are the Floquet 
states momenta and C a (to) are the expansion coefficients of the initial wave function in the 
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FIG. 1: Left panel: (a) Two bands of quasienergies in the linear regime whose interaction leads to 
avoided crossings. Symbols placed in different bands correspond to states with Husimi functions 
depicted in the right panel. Here E2 = 1.2. (b) Average current J vs 6 for different amplitude 
values of the second harmonic E2: 1 (dashed line) and 1.2 (solid line). Right panel: Husimi 
functions of two different Floquet states labeled by a circle and square for to = 0. We use the 
Husimi function defined in (|. Bands and current are plotted for the interval (— 7r,0). Here the 
quasienergies are symmetric with respect to a reflection in 9 as deduced from Eq. tjlip . while from 
Eq. fllOj) it follows that the current is antisymmetric with respect to a reflection at the origin. The 
other parameters are E\ = 3.26 and u = 3. 



basis of Floquet states. Breaking the symmetries S a ([5]) and §6§, we desymmetrize the 
Floquet states, i.e. the Floquet states momenta acquire a finite value {p) a 7^ 0, which results 
in the appearance of a directed transport. 

In general, the current is a function of the initial time to an d the relative phase 9, namely 
J(to, #)• After averaging over the initial time it exhibits the property 



j(6) = -j(9 + n) = -J(-e). 
5 



(10) 



We focus the analysis on previous computations obtained for fi = 0.2 in 




Fig{Tk shows a section of the quasienergy spectrum of the Floquet states versus 9 for the 
interval [— tt,0], where a resonance of states, i.e. an avoided crossing, takes place. The 
quasienergies of the system Eqs. ([2111]) in the linear limit possess the property jl], 0] 

Ue) = ~e a {-6). (11) 

FigHb shows the current dependence on 9. The computation is performed taking the initial 
state |0) = 1/ V2tt, which overlaps with states in the chaotic layer. With this initial condition 
we mimic a dilute gas of atoms which are spread all over the lattice with zero momentum. 
The current has two peaks which are linked to avoided crossings displayed in Figfjji. In these 
particular avoided crossings, states from the chaotic layer and transporting states mix, which 
leads to a leakage from the chaotic layer to the transporting state, thereby enhancing the 
current. 

On the other hand, it was shown in that by tuning the amplitude of the second 

harmonic of the driving force the peaks becomes broader which makes it easier resolving it 
in experiments (see dashed line in FigUJi). 



B. Nonlinear regime 



In the nonlinear case the analysis of the generation of directed transport in the presence of a 



driving 
mixing 



orce is much more complicated, due to possible nonintegrability, classical chaos, and 



151 ] . However, one can take the Floquet states of the linear problem, and continue 



them as periodic orbits into the nonlinear regime. Then these nonlinear Floquet states can 
be analyzed. 

While in the linear regime the evolution of the Floquet state is determined by Eq.(JH]) with 
a unitary operator U, in the nonlinear regime the unitarity is lost, and is replaced locally 



by simplectic maps 



n 



Nevertheless, a similar transformation over one period of the ac 



driving can be defined, viz. 



Uip a (0) = exp(-ie a )-0 a (O), 



(12) 



where U is a nonlinear map of the phase space onto itself, defined by integrating a given 
trajectory over one period of the ac driving [171 ] . The solutions ip a constitute generalizations 
of the linear Floquet states ([7]). 



To compute the nonlinear Floquet states, we use a numerical method implemented in com- 



putational studies of periodic orbits 17|. The basics steps are as follows. First, we choose 
a linear Floquet state as an initial seed. Then taking a small value of the nonlinearity 
strength, we compute the new solution using a Newton-Raphson iterative procedure, by 
variying initial seed (see Appendix A for a detailed explanation). The procedure involves 
conservation of the norm and the variation of the quasienergy of the state, which together 
enforce the convergence to the desired solution. In each iteration step, the new trial solution 
is integrated over one time period T. Once a solution is found, we increase the nonlinear- 
ity strength again by a small amount and repeat the same procedure. Thus, we trace the 
solution into the nonlinear domain. 

Desymmetrization of the Floquet states, due to breaking of symmetries, leads to the appear- 
ance of directed transport in the linear regime with enhancement of transport due to resonant 
Floquet states. It is therefore worthwhile to investigate what happens with nonlinear Flo- 
quet states in the absence of those symmetries, and trace the fate of the abovementioned 
resonances in the nonlinear regime. 

1. Dimer 

To gain insight into the effect of breaking symmetries on nonlinear Floquet states, we utilize a 
basic model of two coupled BEC states in the presence of an external driving. The equations 
for a driven two sites model can be written as 

%H ^ = C>fo + g1> 1 N 1 + il> 1 f(t), (13) 
ifi ^ = CV>i+^ 2 iV 2 -^ 2 /(t), (14) 

where fit) = / 1 sin(a;t) + / 2 sin(2u;i + 9), C is the coupling term, and Aq j2 = IV'i^I 2 are 
the populations or number of particles in the sites 1,2. The above equations can be also 
qualitatively viewed as a restriction of the original case (j3J) to just two basis states with 
opposite momenta. That leads to the corresponding different signs of the last terms in the 
rhs of the above equations. 

With f(t) = the equations above are used to describe, on the mean field level, the self- 



trapping transition of two coupled BEC. It was shown in 18|, [l9| that such a phenomenon 



occurs when the nonlinearity exceeds a critical value, and the new states are characterized 



by a population imbalance N± — No. The existence of critical or threshold values have 

n 

been predicted as well |18j . The critical value is determined by the bifurcation point of the 
stationary solutions. The selftrapped states violate the permutational symmetry of (TT51) . (TT3j) . 
which implies that the equations are invariant under permutation of the two indices. If 
f(t) 7^ 0, the permutational symmetry is broken in general. If however f(t) is symmetric, 
then the dimer equations are invariant under the combined action of permutation and time 
reversal. 

Let us consider the linear case g = first. For f(t) = the stationary solutions are the 
in- and out-of-phase modes ipi(t) — i^CO = ^75 e=R ~*! which are in fact strictly time- 
periodic states, with period 27r/z/C. Adding the ac drive f(t) they transform into two linear 
Floquet states. They will start to become asymmetric, since the presence of f(t) violates 
permutational invariance, but only if f(t) is also violating time-reversal antisymmetry. Note 
that if f(t) is antisymmetric, the original ac field E(t) is symmetric. If on the contrary time- 
reversal antisymmetry is in place, the linear Floquet states will still be invariant under the 
combined action of permutation, time reversal, and complex conjugation. Thus they do not 
acquire any population imbalance. Also the Floquet states will now acquire a nonzero phase 
shift e~ luT ', where v « — in the limit f{t) — > 0, when iterated over one period of the driving 
T. 

In the presence of nonlinearity close to these linear Floquet states there will be again states, 
which are 'periodic' in the sense that after one period of the ac driving the state returns to 
itself, up to a corresponding phase. Note that both the shape of the eigenstate, but also 
the phase (i.e. the analogue of the quasienergy in the linear case) will smoothly change 
upon continuation into the nonlinear regime. The continuation process of such a nonlinear 
Floquet state is thus encoded by the linear Floquet state at g — 0, which is chosen to be 
continued. 

In the case of fi, fi being nonzero, 9 becomes a relevant parameter. For 9 7^ 0,tt time- 
reversal antisymmetry is broken, and the nonlinear Floquet solutions will also loose that 
symmetry, together with the permutational symmetry (see above). 

Hereafter, we name by symmetric the case when 9 = 0, tt and nonsymmetric otherwise. We 
will also coin linear and nonlinear Floquet states periodic orbits, although they are only 
periodic up to the above mentioned phase shift. 

Fig. [2] shows the result of continuation of one periodic orbit for the symmetric and non- 
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FIG. 2: Nonlinearity parameter g vs quasienergies for one of the nonlinear Floquet states, (a): 
6 = 0. (b): 9 = —1.6. Right-down inset shows the average populations over a period T as a 
function of the nonlinearity strength g. Left-up inset is the enlargement around g = of the 
right-down inset. The parameters are u = 2ir, fi = f2 = 1, C = 1, fj, = 1. 



symmetric driven dimer, which corresponds to the symmetric eigenstate of the linear un- 
driven system. The symmetric case shows a pitchfork bifurcation with the appearance of 



two new states at the bifurcation point (Figl2K • 



2d, 
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'hese new states exhibit a population 



22] , despite the fact that the equations 



imbalance similar to the undriven system |18|, 
are symmetric. 

This bifurcation is a hallmark of the presence of nonlinearity. Out of a given Floquet state 
several new states are emerging, while for the linear case the number of linear Floquet states 
is fixed by the size of the chosen basis. 

Conversely, if the time-reversal symmetry is broken, a saddle-node bifurcation appears. First 
of all the state continued from the linear limit, already acquires some nonzero population 
imbalance, since the symmetry is broken. Figj2] shows that the strict continuation of that 
state evolves into a state with a strong population imbalance. Two other states - one which 
is corresponding to a weak imbalance, and another which has a strong imbalance as well - 
emerge through the saddle-node bifurcation. 

To conclude this part, we may expect that nonlinearity induces Floquet states with nonzero 
population imbalance via bifurcations, or enhances the already present imbalance (originat- 
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ing from a symmetry breaking) again via bifurcations. We remind the reader here, that 
the simple dimer model can be obtained from by choosing two momentum basis states 
with different sign, and replacing the original complicated interaction which is mediated via 
further basis states, by a direct interaction term. A population imbalance thus means a 
momentum imbalance as well, i.e. a nonzero current. It may be thus as well expected, that 
for the full problem, to be treated below, nonlinearity may enhance directed currents via 
bifurcations. 

2. Full lattice 

The analysis of nonlinear Floquet states in the dimer suggests the formation of states with 
nonzero mean momentum above a certain threshold value. The question then arises, to what 
extent those features, exhibited by the Floquet states in the dimer, are manifested in the 
full lattice Eqs.pTJ). 

On the other side, the question about the fate of the resonances found in the linear limit 
remains. It is therefore of particular interest to analyse those resonant Floquet states, which 
lead to an enhancement of transport in the linear regime. By switching on the nonlinearity 
in Eq.p]), the system in the linear regime gets perturbed and the resonances are shifted to 
different values of the control parameters. 

In FigJSl we present the continuations of two Floquet states depicted in FigHJ into the 
nonlinear domain. We observe a bifurcation of one of the periodic states as we increase the 
nonlinearity strength. Here, as in the dimer, a bifurcation of saddle-node type leads to the 
formation of three periodic solutions out of one. 

Since we want to study the transport of a BEC with zero initial momentum, we focus on the 
continuation of the state lying in the chaotic layer. The continuation process is indicated by 
a sequence of Husimi functions in FigJHh, and can be summarized as follows: i) before the 
bifurcation, the state is located in the chaotic layer; ii) after the bifurcation, the periodic 
state transforms into a mixed state due to a resonance with a second transporting state; 
iii) further increase of the nonlinearity strength transforms it into a transporting state. By 
contrast, the originally transporting state does not experience significant changes as we 
increase the nonlinearity strength. 

A quantitative measure of this transition process is given by the evolution of the average 
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FIG. 3: (a): Nonlinearity strength vs quasienergy of periodic solutions. The periodic solutions are 
a continuation from the Floquet states displayed in FigfTJ The symbols link to states displayed 
in FigfTJ Insets: Husimi functions {P vs X) (cf. Fig(T]). The sequence of Husimi functions for 
the red line is explained in the text, (b): Enlargement of the region where the quasienergies of 
the continued original states intersect in panel (a). The results from the computation of Eq. (|16p 
appear superimposed to the circles with solid lines. Inset: Mean average momentum of periodic 
solutions vs nonlinearity strength, blue square: momentum of the state depicted by a blue circles- 
dashed line in (a), red circle: momentum of the state depicted by a red circles-dashed line in (a). 
The parameters are E\ = 3.26, E\ = 1.2, 9 = —1.6 and to = 3. 

momentum of the periodic states as we change the nonlinearity strength [FigJ2fc>]- The 
momentum of the state initially located in the chaotic layer goes from small to high values 
with a sharp increase at g ~ 0.003 (i.e. at the bifurcation point), whereas the momentum 
for the transporting state remains nearly a constant. 

So far, we have analyzed the evolution of two states in the nonlinear domain. A clear message 
from the above results is, that nonlinear Floquet states may drastically change their average 
momentum at bifurcations, which are also visible through sharp changes in the dependence 
of their quasienergies on the nonlinearity parameter. 

In order to estimate the shift of the bifurcation point, and to finally predict the new possible 
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resonance positions with increasing nonlinearity, we use a perturbation approach to estimate 
the dependence e = P(e,g), related to the Aharanov-Anandan phase [23'] (see also 16]). 
Assume that nonlinear periodic solutions and Floquet states are related through a dynamical 
phase X(t), i. e., </>(t) = exp[-iX(t)]4>(t) with A(T) - A(0) = 2kn, k = 0, ±1, ±2, ... 
Then for k = we find (see Appendix B for details) 

e{0,g) = ~e{6) + - C WlHV), (15) 
A 1 Jo 

where (...) = J Q 27r ...dx is the inner product defined in the Hilbert space of Floquet states. 
From Eq. (fT5l) it follows that states with strong nonlinear interaction have large quasienergy 
variations. That allows to resonantly couple states located deep in the well, which exhibit 
strong localization, with other states. These resonant couplings are beyond the dimer model 
discussed above. 

Bifurcations of new states correspond to a coalescence of different families of states, leading 
to 4>j{t) = exp[—i\j(t)](a,j(j)i(i) + bjfa(t)), where fa(t) and fait) are the corresponding 
original Floquet states. The evolution of the quasienergies for the nonlinear periodic states 
is thus given by 

efag) = \ aj \%(6) + fe| 2 e 2 (0) + 9 - f dt{fa\ \fa\ 2 \fa), (16) 

^ Jo 

where a,j and bj are the corresponding weightings of the linear states in the nonlinear state 
expansion (see Appendix B for details). 

In Figj3b, we plot the quasienergy values computed with Eq. flTB]) . using the nonlinear states 
continued from the original linear states depicted in Fig{TJ We find an excellent agreement 
with the full numerical results. 

For weak nonlinearity, the quasienergies of the states depend linearly on g. A simple per- 
turbation expansion allows to take the wavefunctions of the original linear states and its 
respective quasienergies with Eq.( jT5i) . Then from the quasienergy intersection of the two 
states, €ii9,g) = €2(0, g), we compute the critical value of g: 

Iq T Jo 7 " \ fa\*dxdt - J'q J Q 27r \fa\ 4 dxdt 

Inserting the wave function of the original linear states depicted in Fig{T]in Eq.l fPTI) . we obtain 
g = 0.003008, which is a good estimate of the nonlinearity strength at the bifurcation point. 
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FIG. 4: Left panel: Current dependence upon the nonlinearity strength g for 8 = —1.2, —1.1, —0.99 
with initial time to = 0. Right panel: (a) Current dependence upon 9 for g=0.005 with to = 0. 
(b) Current as a function of the initial time for 6 = —1.2 and two different nonlinearity strengths. 
Circles: g=0.005, squares: g=0.001. The other parameters are E\ = 3.26, E\ = 1.2, u = 3. 

Let us now investigate the evolution of the state |0) for k — 0. In the linear regime, 

we use the Floquet representation to derive the expression for the asymptotic current. In 

the nonlinear regime it is no longer possible. Instead, we compute the running average 

rt 

momentum P = ^ , \ / p dt with p = (tp\p\tp), over long times, which becomes the 
I* _ *oJ J t0 

asymptotic current in the limit t — > oo, i.e., J (to) = lim^oo P. To validate the convergence 
to the asymptotic current we use maximum integration times ranging from 100000 till 300000 
time periods. 

First we compute the current for 9 pa —0.99. In such a case, for g = 0, there is an avoided 
crossing (see above discussion), i.e. a resonance. Further increase of the nonlinearity leads to 
a decay of the current (FigHj). However, taking 9 = —1.1, —1.2, we observe a corresponding 
shift of the resonance peak to larger values of the nonlinearity parameter, thus showing a 
robustness of the resonant current enhancement observed in the linear limit. 
FigJH shows the dependence of the current on 9. Resonances similar to the linear case are 
observed, which are shifted in 9. Thus the nonlinearity strength can also be used as a control 
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parameter for tuning resonances. 

It was shown for the linear regime that the current depends upon the initial time 

tn- In the nonlinear case, the dynamics may exhibit similar behavior to classical chaos 



15j, making the analysis more complicated. An essential point here is that, while in the 



linear regime Floquet states mix in narrow parameter regions via avoided crossings, in the 
nonlinear domain mixed states may survive for a fairly large range of nonlinearity values. 
This, along with the fact of having multiple bifurcations with the ap pea rance of new periodic 
states [241 ]. may lead to classical-like chaotic behavior predicted in [15 ]. 



A benchmark for the persistence of directed transport is that the sum of currents for different 
initial times do not cancel. To check that, we compute the momentum evolution for different 
initial times with the system in and out of resonances. In FigJH we present the current 
dependence upon the initial time to. The curve depicted by filled circles shows a similar 
behavior to the one obtained in the linear regime with the system in resonance (cf. Fig. 4 
in [5I). It also displays a large positive current for all initial times, thus confirming the 
existence of directed motion. On the contrary, the curve with squares shows a nearly flat 
profile with smaller current values, recalling the off resonance scenario in the linear limit. 



III. CONCLUSIONS 



We have studied the rectification of interacting quantum particles in a periodic potential 
exposed to the action of an external ac driving, using a mean-field approach. 
We showed that by tuning the nonlinearity in an optical lattice it is possible to enhance the 
directed transport of cold atoms. A possible experimental way to achieve it is to vary the 
scattering length by changing the strength of the magnetic field [251 ] . These resonances are 
partly continued from the noninteracting system, but become saddle node bifurcations of 
more complicated nonlinear Floquet states in the nonlinear system under consideration. 
We developed analytical estimates of the nonlinearity strength in the resonance, and showed 
that the evolution of an initial state with zero momentum carries all signatures of a ratchet 
state, i.e. its average momentum is nonzero. Therefore, ratchets and quantum resonances of 
(nonlinear) Floquet states are robust with respect to interactions, and should be observable 
in real experiments. 

Finally, bifurcations of the nonlinear Floquet states have been observed, which may affect 
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the measured currents strongly. 



APPENDIX A: NUMERICAL PROCEDURE FOR THE COMPUTATION OF NON- 
LINEAR FLOQUET STATES 

Take Eq.(fT2D 

Uif) a (0) = exp(-ie a )V>«(0), (Al) 

where ip a are periodic solutions in the nonlinear domain and U is a symplec- 
tic operator. The periodic states ^ Q (0) are decomposed into real and imagi- 
nary parts as = {Real[ip a (0)}, lm[ip a (Q)]}. Likewise, we define the vector Q = 
{Real[exp(ie a )U-ip a (0)},lm[exp(ie a )U'ip a (0)}}. For convenience we write them as 2N dimen- 
sional vectors X a (0) = $ Q (0), X a = ^[^«(0)], where N is the dimension of the Hilbert space 



261 ] . Hereafter, for simplicity, we drop the index a. 

The two vectors X, X(0) are identical for the linear case. Assuming weak nonlinearity, we 
compute the vector X after one integration period, taking as initial seed the linear state 
X(0). This implies that after a full integration period our final state deviates from the initial 
state: 

G[X(0)}=X -X{0). (A2) 

To correct such a deviation we use the Newton-Raphson method. Basically, the method 
updates the initial seed X(0) after every iteration until each component of G is reduced to 
a value less than 10~ 9 . For every integration over the period T, we use the split operator 
method. 

To successfully accomplish the reduction of the vector difference ( 1A2I) . some constraints 
should be fulfilled. First, the Floquet states should preserve the norm. We thus add an- 
other component G2N+1 = XX — X(0)X(0), and the iteration process is also zeroing this 
component. 

To fully characterize our nonlinear Floquet states, we define the new variable vector Y = 
(X,e), with the quasienergy as the 2N + 1 component. Next, we define the new vector 
function F(Y) = [G(Y),G2N+i{Y)]. Having 2N + 1 variables and 2N + 1 functions to be 
zeroed, we now may apply standard Newton-Raphson methods. 
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APPENDIX B: QUASIENERGY DEPENDENCE ON THE NONLINEARITY 
STRENGTH 



We define an expression similar to Eq.ffl), but for the nonlinear case, viz. 



\m) = e-**Wf)>, m+T)) = \m- ( Bi ) 



Here, we have dropped the index a. 

Then, inserting Eq. flBll) into Eq.([2]), one finds 



i^\<t>)=H\<t>), (B2) 



where Ti = H — + g\(fi\ 2 . Similarly, using Eq. ([7]) we get for the linear case 

i^\$)=H\4>), (B3) 

where H = H — 

Now define <f>(t) = exp[-i\(t)]4>(t) such that A(T) - A(0) = 2kir with k = 0, ±1, ±2... This 
assumes that after projecting the nonlinear state onto the linear space, both vectors, the 
projected and the Floquet states in the linear space, have a phase difference X(t) [3]. For 
simplicity we take in the following k = 0. 
Using then <j>(t) and Eq.( lB3l) . we get 

Using Eq. (1B2I) and integrating over a period, we find 

1 ' T 



= - / dt[(<P(t)\H Q \<j>(t)) - (<!>(t)\H \<j>(t))] 
V Jo 



-~e + e- 9 - [ T dt(m 2 W- (B5) 
V Jo 



The first term of the r.h.s. of equation (1B5I) is zero after integrating over one period. Then 
we obtain 



e = ~e+ 9 - [ T dt(m 2 \&- (B6) 
V> Jo 

As follows from Eq. (1A1|) . e denotes the phase accumulated after a completion of one period 
T . To understand the physical meaning of Eq. (1B6|) . we rewrite it as a function of the real 
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quasienergy £ (see [27] for definition). It relates to our e as £ = ifi/T, where —fx u>/2 < 
£ < \i uj/2. Then we obtain 

£ = £ + ^ J* dt(<P\\<P\ 2 \<P) = £ + ^ J\t gj^dx |0| 4 . (B7) 

This expression tells us that the energy of our nonlinear Floquet state is the sum of the 
corresponding linear and nonlinear contributions. The second term is the time average of 
the state energy due to the nonlinear interaction. 

So far, we have considered the energy evolution of single states without perturbation. How- 
ever, single states may "coalesce" leading to bifurcations. We take the simplest case with 
two resonant states. We project vectors on a basis that results from the superposition of 
the continued original eigenstates, whose continuation lead to resonances in the nonlinear 
domain. That is, 4>j(t) = exp[— i\j(t)}(aj(pi(t) +6j0 2 (t)) with j = 1, 2; where 4>i(t) and fait) 
are the original Floquet states. 
By doing similar operations as above we get 

e s = \ aj \% + \h 3 \ 2 £ 2 + % [ ^11^1%). (B8) 

1 Jo 

where |aj| 2 = |(0i|0j)| 2 , and \bj\ 2 = 1 — |a,-| 2 = |(0 2 |0j)| 2 - 
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